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Abstract. The renewed interest in investigating quaternionic quantum mechanics, in 
particular tunneling effects [1], and the recent results on quaternionic differential opera- 
tors [2] motivate the study of resolution methods for quaternionic differential equations. 
In this paper, by using the real matrix representation of left /right acting quaternionic op- 
erators, we prove existence and uniqueness for quaternionic initial value problems, discuss 
the reduction of order for quaternionic homogeneous differential equations and extend to 
the non-commutative case the method of variation of parameters. We also show that the 
standard Wronskian cannot uniquely be extended to the quaternionic case. Neverthe- 
less, the absolute value of the complex Wronskian admits a non-commutative extension 
for quaternionic functions of one real variable. Linear dependence and independence of 
solutions of homogeneous (right) H-linear differential equations is then related to this 
new functional. Our discussion is, for simplicity, presented for quaternionic second order 
differential equations. This involves no loss of generality. Definitions and results can be 
readily extended to the n-order case. 

PACS. 02.10.Hh - 02.30. Tb MSC. 11R52 - 30G35 - 46S20 



I. INTRODUCTION 

Let K, C = span {1, i}, and HI = span {1, i,j, k} be the real, complex , and quaternionic field [3], 

i 2 = j 2 = k 2 — ijk = — 1 , 

and 

T : R -> R 

be the set of real functions of real variable. Through the paper, quaternionic functions of real variable, 
^>(x) £ H ® J 7 , will be denoted by Greek letter and constant quaternionic coefficients by Roman letter. 
To shorten notation the prime and double prime in the quaternionic functions shall respectively 
indicate the first and second derivative of quaternionic functions with respect to the real variable x, 

df „ d 2 * 

:= — and := - — . 

da; da; 2 

Due to the non-commutative nature of quaternions, it is convenient to distinguish between the left 
and right action of the quaternionic imaginary units i, j, and k by introducing the operators L q and 
R p whose action on quaternionic functions 'J' is given by 

L q ^> = q * and i? p * = *p , (1) 



2 



Stefano De Leo, Gisele C. Ducati: QUATERNIONIC DIFFERENTIAL EQUATIONS 



These (left/right acting) quaternionic operators satisfy 

L q L p = L qpl R q R p = R pq and [L q ,R p }=0, 

and admit for 

q = q + iq 1 +jq 2 + kq 3 , p = p + i p 1 + j p 2 + k p 3 , * = * + i *j + j * 2 + k * 3 , 
the following real matrix representation [4-6] 

/ p„ -p x -p 2 -p 3 ' 



(2) 



( q a -qi -q 2 -q 3 
qi q a -q 3 q 2 
q 2 q 3 q -qi 

\ q 3 -q 2 qi q a , 



Pi Pa P 3 ~P 2 
P 2 -Pa Pa Pi 
\P 3 P 2 -Pi Pa, 



) T . 



(3) 



II. EXISTENCE AND UNIQUENESS. 

In this section we discuss existence and uniqueness for the quaternionic initial value problem 

V'^a^ + p^ + p, *(*„) = /, *'(x )=g, (4) 
with a(x), P(x), p(x) e H ® T, x a e I : (x_,x + ) and /, g e H. 

Theorem 1 . Let a, /3 and p in Eq.(4) be continuous functions of x on an open interval I containing the 
point x — x . Then, the initial value problem (4) has a solution \& on this interval and this solution 
is unique. 

Proof. By using the real matrix representation (3), we can immediately rewrite the quaternionic initial 
value problem (4) in the following vector form 



*2 



/ a 




-Ci 2 


-a 3 ^ 




r* ] 


a x 


Cta 




ct 2 






a 2 


a 3 


Cta 


-a ± 






\0i 3 


-a 2 


a 1 


Cta) 







/ft -A -&-/V 

Pi Pa~P 3 P 2 
P 2 P 3 Pa -Pi 
\P 3 -P 2 Pi Pa, 







' Pa' 






Pi 


*2 


+ 


Pi 


.*s. 




-P 3 . 



(5) 



with 







r/oi 




r*o(*c)i 




'9a' 


*i(ar ) 

y 2 (Xo) 




fi 

u 


and 


*i(*o) 
y 2 (x ) 




9i 

92 


.y 3 (Xa)_ 




./a. 








.g 3 . 



(6) 



Eq.(5) represents a (nonhomogeneous) linear system with a m ,/3 m ,p m £ where m = 0,1,2,3. 

These functions are (see hypothesis of Theorem 1) continuous (real) functions of x on an open interval 
I containing the point x — x a . Then, by a well-known theorem of analysis, see for example ref. [7], 
the linear system (5) has a solution 



on this interval satisfying (6), and this solution is unique 



III. LINEAR INDEPENDENCE AND DEPENDENCE OF SOLUTIONS 

Let us now analyse the linear independence and dependence of the solutions of second order homoge- 
neous differential equations 

*" = a + , (7) 
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where a and j3 are (quaternionic) continuous functions of x on an open interval I. Eq.(7) is linear 
over H from the right. Consequently, if ip is a solution of Eq.(7) only the function obtained by right 
multiplication by constant quaternionic coefficients, ipu, still represent a solution of such an equation. 
The general solution of Eq.(7) is given in terms of a pair of linearly independent solutions tp and £ by 

* = ipu + ^v , (8) 

where ip — (p„ + i ip t + j <p 2 + k ip 3 , £ = £ + i £ x + j £ 2 + k £ 3 € H <g> T and u, v e H. 

In the standard complex theory (ip = <p> a + i <p 1 and £ = £ + i £ x e C <g> JF) a useful criterion to es- 
tablish linear independence and dependence of two solutions of homogeneous second order differential 
equation, uses the concept of Wronskian of these solutions defined by 

W = p£' -<p' £, , WeC®f . (9) 

This definition cannot be extended to quaternionic functions. Let us consider two linearly dependent 
solutions of Eq.(7), 

£ = <pq, ^,(ei®f , qeU. (10) 
By substituting £ = ipq and £' = ip' q in the Wronskian (9), we find 

<p? - = <p<p'q- <p'vq ¥= . 

Observe that a quaternionic function and its first derivative do not, in general, commute. Thus, the 
definition (9), and all its possible factor combinations cannot be extended to the quaternionic case. 

Let us now use the linear dependence condition (10) to investigate the possibility to define a 
quaternionic functional which extends (in a nontrivial way) the standard (complex) Wronskian to the 
noncommutative case. From Eq.(10) and its derivative, we get 

where tp^ 1 = 1/ip and (<p>') _1 = l/<p'. Consequently, for linearly dependent quaternionic solutions , we 
have 

c-^v^o. (ii) 

To recover, in the complex limit, the standard definition (9) we multiply £' — ip'p^ 1 ^ by ip. Due to the 
noncommutative nature of quaternions, we have to consider the following possibilities 

W L = <p(?-tf/<p- 1 Z) and W R = (e-<p'<p- 1 Z) cp. (12) 

Obviously two other similar definitions can be obtained by ip <-» £, 

W L = -Z{tft -Z'S- 1 tp) = -W I \<p~f\ and W R = -{tff-Z'f 1 tp)S = -W R \tp~Z]. (13) 

The quaternionic functional (12) and (13), which give in the complex limit the standard definition, 
extend a first important property of Wronskian. Two solutions of Eq.(7) are linearly dependent on 

I If Hl(s) W l ( r ) is zero on I. To avoid ambiguity in defining the Wronskian, we shall introduce a 

(real) functional, 

\w\> = \w L \> = \w R \* = \w L \* = \w R \* , 

which extends the squared absolute value of the Wronskian. This unique functional is 

iwi a Hd 2 ief +Kiviw p. (w) 

where ip c = tp — i <p 1 — j p 2 — k <p 3 and £ c = £ — i £ x — j £ 2 — k £ 3 are, respectively, the quaternionic 
conjugate functions of ip and £. 

Observe that Eq.(14) can also be obtained as an application of the Dieudonne theory of quater- 
nionic determinants [8-12]. In fact, 

\W\ 2 = [Det(M)] 2 :=det (MM+) , (15) 
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where 



M 



Theorem 2 . Let a and (3 in Eq.(7) be continuous functions of x on an open interval I : (a, b). Then, 
two solutions ip and £ of Eq.(7) on I are linearly dependent on I if and only if the absolute value of 
the Wronskian, \W\, is zero at some x in I. 

The proof will be divided into three steps: 

(a) - If ip and £ are linearly dependent on I then \W\ = 0; 

(b) - If \W\ = at some x in I then \W\ = on I. 

(c) - If \W\ = at some x in I then if and £ are linearly dependent on I. 

Proof (a). If if and £ are linearly dependent on I, then Eq.(10) holds on I. From Eq.(10), we get 

\w\> = M 2 \ipr \q\ 2 + M 2 kT M 2 - M 2 v' c <p' M 2 *>W C = o , 

then |W| = 0. 

Proof (b). Let us consider Eq.(14). By calculating the first derivative of the left and right side term, 
we obtain 

2 \W\ \W\' = <f'ip£C + w?£'Z + W>c*'& + Wc?'e+ 

- *;v.«: - vv:^ - vv.r£ - ^Vc£e+ 

= \f\ 2 (*x + eg) + lerwv: + vv/) + 

- *;v.«: - vVc^e - *^cw: - £'£W C ' 

= M 2 (a |£T + + h.c.) + |£| 2 (a \ipT + /W. + h.c.) + 

- [(a ^Vc + fi M 2 ) it + h.c] - [(a £'£ c + (3 |£| 2 ) ^ + h.c] 
= 2Re[a] (M 2 1£'| 2 + |£| 2 1^'| 2 - <^Vc££ - M) 

= 2Re[a]|iy| 2 . 

By a simple integration, we find 

|W(a;)| =cxp 

This prove the statement (b). 
Proof (c). From the statement (b), we have 

\W{x )\=Q 

This implies that the quaternionic matrix 



Re[a{y)] dy 



\W(x )\ 



(16) 



\W(x)\=o, xei 



is not invertiblc on I [12]. Hence the linear system 

<P Qi + £ ?2 = , 
<P' Qi + €' I2 = , 

in the unknowns q li2 € H, has a solution (q 1 ,q 2 ) where q 1 and q 2 are not both zero. Recalling that ip 
and £ are linearly independent on an interval I if 



<p(x) qi + £0) q 2 = o 



91=92=0 
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the fact that q x and q 2 are not both zero guarantees the linear dependence of ip and £ on I ■ 
Example 1. Show that <p — exp[— ix] and £ = exp[(i — j)x] form a basis of solutions of 

*"+j* / + (l-fc)* = 0, (17) 

on any interval. 

Solution. Substitution shows that they are solutions, 

[-1 + + 1 - k] exp[-ia;] = , 
[-2 j) + l-k] exp[(z - j)x] = , 

and linear independence follows from Theorem 2, since 

\W\ = ^f\i-i\ 2 + \i\ 2 + i{]-i)-{i-])i = V5 . 

IV. HOMOGENEOUS EQUATIONS: REDUCTION OF ORDER 

Let ip be solution of Eq.(7) on some interval I. Looking for a solution in the form 

S,=pT 

and substituting £ and its derivatives 

4 = ip t + ip t and £ = <p t + Zip t + tpr 

into Eq.(7), we obtain 

T " = (p- 1 aip-2p- 1 p') t' . (18) 

It is important to observe that, for quatcrnionic functions, we cannot give a formal solution of the 
previous equation. Only in particular cases, Eq.(18) can be immediately integrated. For example, for 
homogeneous second order equations with constant coefficients, 

a{x) -> a e H and 0(x) -» b e H , 

at least one solution is in the form of a quaternionic exponential, p = exp[gx], and consequently 
Eq.(18) reduces to 

pr" = (a- 2q) tpr' . (19) 
Let us introduce the quaternionic function 

a = tpr' . 

Observing that 

a' = p'r' + pr" = q pr' + pr" , 
Eq.(19) can be rewritten as follows 

a' = (a-q)a. (20) 
This equation can be immediately integrated, its solution reads 

a = cxp[(a — q)x] . 

Thus, the second solution of the homogeneous second order differential equation with constant coef- 
ficients is given by 

£ = exp[qx] J exp[— qx] exp[(a — q)x] Ax . (21) 

In the complex limit (a, q € C) we find the well-known results £ oc cxp[(a — q)x] if 2q ^ a and 
£ oc xexp[ga:] if 2q = a. In the quatcrnionic case (a,q e C), the integral which appears in (21) 
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must be treated with care. The solution of this integral will give interesting information about the 
second solution of quaternionic differential equations with constant coefficients when the associated 
characteristic quadratic equation has a unique solution. To solve the integral in Eq.(21), we start by 
observing that 

[e ux e vx }' = ue ux e vx + e ux e vx v = (L u + R v ) e ux e vx . 
If the operator L u + R v is invertible the previous equality implies 



/ 



e ux e vx dx= (L^ + RX 1 e ux e vx 



This result guarantees that, if the operator L_ q — R a _ q is invertible the second solution can be written 
in the form 

£ = cxp[gx] {L_ q + i? a _J _1 exp[-gx] exp[(a - q)x] 
= cxp[qx] {R a _ q - L^ 1 cxp[-<?x] exp[(a - q)x] . (22) 

If the operator L_ q + i? a _, is not invertible, we need to solve the integral which appears in (21) by 
using the polar decomposition of quaternions (see example 3) and a term linearly dependent on x 
will appear. In the complex case (a,q G C), the operator L_ q + R a _ q is not invertible if and only if 
2q = a. In the quaternionic (a, ggl), the condition 2q ^ a does not guarantee that the operator is 
invertible. 

Example 2. Knowing that (p = cxp[— ix] is solution of the homogeneous second order equation (17), 
find (by using the method of reduction of order) a second independent solution, £. 

Solution. We have q ~ —i and a = —j. To use Eq.(22) we have to prove that the operator 

L_ q + R a _ q = Li + Ri-j 
is invertible. A simple algebraic calculation shows that 

(L 4 - {L t + // , = 1 . 

Thus, 

{L^ + R^)- 1 =L i -R i _ j . 
We are now ready to calculate £ from Eq.(22), 

£ = cxp[— ix] {L t — Ri-j) exp[ix] exp[(i — j)x] = {Li — Ri-j) exp[(i — j)x\ = exp[(i — j)x] j . 

Due to the H linearity (from the right) of Eq.(17) the right factor j can be ignored recovering the 
solution of example 1. 

Example 3. Inspection shows that 

• / 'I'' • 4 I) (23) 
has <p = exp[— x] as a first solution. Find the second linear independent solution. 

Solution. We have q = — and a = —i. In this case, the operator 

L_ q + R a _ q = Ltp, + R^i 

is not invertible. This is easily seen by using, for example, the real matrix representation (3). Thus, the 
integral in Eq.(21) cannot be expressed in terms of an exponential product. Let us explicitly calculate 
£ from Eq. (21). Wc find 

i = cxp[-^2- x] J exp[^ x] exp[^ x] dx 

= exp[-^z] J (cos^ + ^sin^) (cos^ + ^sin^) dx 

= cxp[-^x] J {1 - k exp[-(i + j)x]} ^ dx . 
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Due to the H linearity (from the right) of Eq.(23) the right factor can be removed. After integra- 
tion, we find 

£ = cxp{-i±2-x] {x- k 1 -^- exp[-(i + j)x}} = (x + cxp[-^x] . 

Observe that the quaternionic factor appears on the left of the quaternionic exponential and 
consequently cannot be removed. It is a fundamental part of the solution. Inspection shows that 

£ = xcxp[—^-x] 

is not solution of Eq.(23). 



V. NONHOMOGENEOUS EQUATIONS: VARIATION OF PARAMETERS 

A general solution of the nonhomogeneous equation (4) is a solution of the form 

* = * h + * P , (24) 

where 

4\ = tp q x + £ q 2 

is a general solution of the homogeneous equation (7) and ^ p is any particular solution of (4) containing 
no arbitrary constants. In this section we discuss the so-called method of variation of parameters to 
find a particular solution for quaternionic nonhomogeneous differential equations. 

A method to solve a homogeneous second order quaternionic differential equations with constant 
coefficients has been recently developed [2]. Quaternionic differential equations with non constant 
coefficients are under investigation. We suppose to know two independent solutions of the homogeneous 
equation associated with Eq.(7). We wish to investigate if the method of variation of parameters still 
works in the quaternionic case. 

The method of variation of parameters involves replacing the constant q 1 and q 2 by quaternionic 
functions i' 1 (x) and v 2 (x) to be determined so that the resulting function 

is a particular solution of Eq.(4). By differentiating * p we obtain 

The requirement that ^f p satisfies Eq. (4) imposes only one condition on v x and u 2 . Hence, we can 
impose a second arbitrary condition, that is 

^I+£^=0. (25) 

This reduces <£' to the form 

V 

By differentiating this function we have 

Substituting * p , <J p , and in Eq.(4) we readily obtain 

^X+CX = p. (26) 
Collecting Eq.(25) and Eq.(26), we can construct the following matrix system 

(27) 









"o" 








p 
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from which {\W\ ^ 0) we obtain 





-(yl) 


0" 


<_ 




p 



[v-^'-V]" 1 



^-rrv] -1- 



) 


"o" 







Then, 



^ = b'-£'CVrV and i^=[^-^V- 1 €]" 1 p. 
To find ^i(x) and ^ 2 (x) we have to integrate the previous equations. 

Example 4. Find a general solution of the nonhomogeneous quaternionic differential equation 

+ jV + = ix 

Solution. The solution of the associated homogeneous equation (see example 1) is 

* h = exp[-ia;] q t + exp[-(i + j)x] q 2 . 

The particular solution is 

* p = cxp[-za;] v x + exp[-(i + j)x] v 2 . 
Consequently, from Eqs.(29) we find 

v[ = exp[za;] x k and v' 2 = — cxp[(z + j)x] x k 

which after integration give 

i^i (x) = (1 — ix) exp[ix] k and u 2 (x) — —\ [1 — (i + j)x] exp[(i + j)x] k . 

Finally 

*p = 5 [(i+j)x + k] . 

A general solution of Eq.(23) is 

* = cxp[-za;] q 1 + exp[— (i + j)x] q 2 + \ [(i + j) x + k] . 



(28) 
(29) 



(30) 



VI. CONCLUSIONS AND OUTLOOKS 

The recent results on violations of quantum mechanics by quaternionic potentials [1] and the possibility 
to get a better understanding of CP-violation phenomena within a quaternionic formulation of physical 
theories [1,13] stimulated the study of quaternionic differential operators [2]. In this paper, we have 
proved existence and uniqueness for quaternionic initial value problems and solved simple quaternionic 
differential equations by discussing the reduction of order for quaternionic homogeneous equations and 
by extending to the non-commutative case the method of variation of parameters and the definition 
of absolute value of the Wronskian functional. 

In view of a more complete discussion of quantum dynamical systems using quaternionic wave 
packets, our next research (mathematical) interest will be the study of quaternionic integral trans- 
forms. The quaternionic formulation of Fourier transforms could find an immediate and interesting 
application in the study of delay time modifications of wave packets scattered by a quaternionic 
potential step. 



Acknowledgements 

The authors thank the referee for comments, references, and suggestions, and for drawing their atten- 
tion to an alternative way to obtain the squared absolute value of the Wronskian (14) by Diuedonne's 
theory of determinants [see Eq.(15)]. 



Stefano De Leo, Gisele C. Ducati: QUATERNIONIC DIFFERENTIAL EQUATIONS 



9 



References 

1. S. De Leo, G. C. Ducati and C. C. Nishi, Quatcrnionic potentials in non relativistic quantum mechanics, 
J. Phys. A 35, 5411-5426 (2002). 

2. S. De Leo and G. C. Ducati, Quaternionic differential operators, J. Math. Phys. 42, 2236-2265 (2001). 

3. W. R. Hamilton, Elements of Quaternions, vol I, 3 rd Ed, Chelsea Publishing Co., New York (1969). 

4. S. De Leo and G. C. Ducati, Quaternionic groups in physics, Int. J. Theor. Phys. 38, 2197-2220 (1999). 

5. S. De Leo and G. Scolarici, Right eigenvalue equation in quaternionic quantum mechanics, J. Phys. A 33, 
2971-2995 (2000). 

6. S. De Leo, G. Scolarici and L. Solombrino, Quaternionic eigenvalue problem, J. Math. Phys. 43, 5815-5829 
(2002). 

7. E. Kreyszig, Advanced Engineering Mathematics (John Wiley & Sons, New York, 1993), pag. 165. 

8. E. Study, Zur Theorie der linear Gleichugen, Acta Math. 42, 1-62 (1920). 

9. J. Brenner, Applications of the Diuedonne determinant, Linear Algebra Appl. 1, 511-536 (1968). 

10. F. J. Dyson, Quaternion determinants, Helv. Phys. Acta 45, 289-302 (1972). 

11. H. Aslaksen, Quaternionic determinants, Math. Intelligencer 18, 57-65 (1996). 

12. N. Cohen and S. De Leo , The quaternionic determinant, Elec. J. Lin. Alg. 7, 100-111 (2000). 

13. S. Adler, Quaternonic Quantum Mechanics and Quantum Fields (Oxford University Press, New York, 
1995). 



